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Abstract 

A  physics-based  single  particle  model  was  used  to  simulate  the  life  cycling  data  of  a  lithium  ion  cell.  The  simulation  indicates  that  there  are 
probably  three  stages  of  capacity  fade  in  a  lithium  ion  cell  used  at  low  rates.  In  the  first  stage,  lithium  ions  are  lost  to  a  film  formation  reaction  (e.g. 
SEI  formation)  and,  consequently,  the  cathode  becomes  less  intercalated  during  cycling.  In  the  second  stage,  the  loss  of  active  cathode  material 
outpaces  the  loss  of  lithium  ions  and  the  cathode  gradually  becomes  more  intercalated  at  the  end  of  discharge.  The  anode  is  the  limiting  electrode 
in  stages  one  and  two  and  the  change  in  the  anode  voltage  causes  the  cell  to  reach  end  of  discharge  voltage.  In  the  third  stage,  the  limiting  electrode 
shifts  from  the  anode  to  the  cathode,  and  the  anode  becomes  increasingly  less  discharged  at  the  end  of  discharge.  Thus,  more  and  more  “cyclable” 
lithium  ions  are  left  inside  the  anode,  which  causes  additional  capacity  fade. 

©  2008  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Secondary  lithium  ion  batteries  are  being  actively  researched 
for  the  application  into  aerospace,  medical,  military  and  trans¬ 
portation  applications.  Battery  life  prediction  is  an  important 
and  active  research  area.  A  substantial  amount  of  work  [1-11] 
has  gone  into  understanding  the  capacity  fade  and  predicting 
battery  life.  For  example,  Stamps  et  al.  [6]  studied  the  capacity 
fade  of  a  lithium  ion  cell  over  1600  cycles  using  a  hybrid  estima¬ 
tion  algorithm  combining  elements  of  both  batch  estimation  and 
online  moving  horizon  estimation.  However,  the  battery  model 
used  in  their  work  was  a  simplified  empirical  model  which  does 
not  provide  insight  concerning  capacity  fade  mechanisms.  Gang 
et  al.  [7]  used  a  physics  based  model  to  simulate  the  cycling  per¬ 
formance  of  lithium  ion  cells.  The  SEI  formation  on  the  anode 
during  charge  was  assumed  to  account  for  the  capacity  fade. 
Lee  et  al.  [8]  simulated  the  cycling  performance  of  lithium  ion 
cells  under  low  earth  orbit  conditions  using  the  same  model. 
Good  agreement  was  achieved  between  the  experimental  and 
simulated  discharge  capacities.  Other  researchers  [11]  also  used 
equivalent  circuit  models  (ECM)  to  study  the  capacity  fade  of 
lithium  ion  cells. 
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In  order  to  understand  how  the  capacity  fade  in  individual 
electrodes  would  affect  the  capacity  fade  at  the  cell  level,  we 
used  a  physics-based  single  particle  model  [9,10]  to  analyze  the 
cycling  data  of  a  lithium  ion  cell.  In  comparison,  most  empirical 
methods  such  as  experimental  data  extrapolation  view  the  cell 
as  a  whole  entity  and  could  not  provide  us  such  information. 

The  focus  of  this  work  is  to  show  a  significant  findings  based 
on  our  data  analysis  on  the  cycling  data  of  a  lithium  ion  cell,  but 
not  to  develop  a  new  cell  model.  Our  study  reveals  the  capac¬ 
ity  fade  of  the  cell  could  have  three  stages  under  low  rate.  The 
experimental  data  show  the  first  two  stages,  controlled  mostly 
by  loss  of  active  lithium  ions  and  loss  of  active  cathode  mate¬ 
rial,  respectively.  A  third  stage  with  accelerating  capacity  fade, 
as  predicted  by  our  simulation,  would  happen  when  the  limit¬ 
ing  electrode  shifts  from  the  anode  to  the  cathode.  The  detailed 
discussion  on  this  capacity  fade  sequence  can  be  found  in  Sec¬ 
tion  3  of  this  paper.  Although  our  experimental  data  have  not 
yet  confirmed  this  three  stage  capacity  fade  pattern,  Bloom  et 
al.  [4]  observed  a  similar  pattern  in  their  cycle  life  study  (Fig.  3 
in  Ref.  [4]).  However,  without  further  detailed  analysis,  it  is  not 
suggested  that  our  findings  can  be  applied  directly  to  their  study. 

To  our  best  knowledge,  there  has  not  been  any  publication 
(in  terms  of  simulation)  in  the  literature  to  show  a  similar  three- 
stage  capacity  fade  pattern  as  ours.  We  also  provide  a  detailed 
and  plausible  explanation  on  the  reasons  behind  the  formation 
of  the  three-stage  pattern.  In  addition,  our  findings  highlight 
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Nomenclature 

Cq  Li+  concentration  in  the  electrolyte  (mol  cm-3) 
cmax  maximum  Li+  concentration  in  the  particles 
(mol  cm-3) 

cs  Li+  concentration  in  the  solid  particles 
(mol  cm-3) 

cs  volume  averaged  Li+  concentration  inside  spher¬ 
ical  particles  (mol  cm-3) 

Ds  solid  phase  Li+  diffusion  coefficient  in  the  parti¬ 
cles  (cm2  s-1) 

4Pp  applied  current  (A) 

jint  intercalation  current  density  (A  cm-2) 

k  kinetic  rate  constant 

n  number  of  terms  in  Eq.  (5) 

Rs  radius  of  the  spherical  particles  (cm) 

U  open  circuit  potential  of  the  electrode  (V) 

V  volume  of  the  electrode  (cm3) 

v  the  stoichiometric  number  at  the  discharged  state 

xo  the  stoichiometric  number  at  the  charged  state 

Greek  letters 

s  volume  fraction  of  the  active  material  in  the  elec¬ 

trode 

km  positive  eigen-value  determined  from 
km  =  tan(A,m). 

0  electrode  potential  (V) 

Subscripts 

i  positive  or  negative  electrode 

neg  negative  electrode 

pos  positive  electrode 


the  importance  to  consider  the  interaction  between  individual 
electrodes  and  the  cell  in  making  the  cycle  life  predictions.  We 
believe  these  contributions  could  provide  useful  information  for 
battery  researchers  to  consider  in  their  own  studies. 

2.  Model 

The  experimental  discharge  capacities  are  presented  in  Fig.  1 
for  a  lithium  ion  cell  (lAh)  with  a  LiNiCo02  cathode  and  a 
carbon  anode.  The  discharge  current  was  8  mA  (<C/100).  The 
discharge  data  were  studied  using  a  so-called  single  particle 
model  [9,10].  The  model  uses  the  assumption  that  the  behavior 
of  the  porous  electrodes  can  be  represented  by  spherical  particles 
of  the  active  materials  in  the  electrodes.  The  concentration  and 
potential  distributions  in  the  solution  phase  are  assumed  to  be 
negligible,  which  is  a  reasonable  assumption  under  low  rate  situ¬ 
ations.  Therefore,  the  usage  of  the  single  particle  model  is  limited 
to  low  current  situations.  The  detailed  description  and  discus¬ 
sion  of  the  single  particle  model  can  be  found  in  the  literature 
[9,10]. 

The  model  equations  consist  of  the  solid  phase  diffusion 
equation  (Fick’s  diffusion  law)  in  a  spherical  particle  and  the 


Fig.  1.  The  experimental  discharge  capacities  vs.  cycles. 
Butler  Volmer  equation  for  both  electrodes: 
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The  intercalation  current  density  j int,/  is  calculated  based  on 
the  equivalent  electro-active  surface  area  (AO  of  the  electrode 
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where  V;  is  the  volume  of  the  electrode  and  £[  is  the  volume 
fraction  of  the  active  electrode  material.  The  decrease  of  the 
volume  fraction  St  could  indicate  a  loss  of  active  material  in  the 
electrode  (i.e.,  through  isolation  of  active  material  particles  from 
the  electrode). 

The  diffusion  equation  (Eq.  (1))  can  be  further  simplified 
from  a  partial  differential  equation  to  a  differential  and  algebraic 
equation  using  the  pseudo  steady  state  (PSS)  method  [12,13]. 
The  PSS  equations  for  the  diffusion  in  the  electrode  are: 
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where  the  number  of  terms  “n  ”  in  the  summation  is  typically  set 
equal  to  5;  csj  the  volume  averaged  concentration  in  the  spher¬ 
ical  particle,  and  csj\r=Rsi  is  the  concentration  at  the  surface 
of  the  particle.  The  state  of  charge  (SOC)  or  the  stoichiometric 
number  of  the  electrode  is  defined  as 


cmax,z 

The  km  values  are  determined  from  the  eigen- value  equation: 
Xm  =  tan  (km)  (7) 

3.  Results  and  discussion 

3.1.  Parameter  estimation 

Three  model  parameters  were  found  to  change  with  cell 
capacity  fade  when  fitting  the  discharge  data.  They  are  the  initial 
SOCs  for  the  cathode  and  the  anode  at  the  beginning  of  charge 
(T0,pos  and  vo,neg)>  an(*  the  volume  fraction  of  the  cathode  (£pos). 
The  volume  fraction  of  the  anode  does  not  change  much  with 
cycling.  Fig.  2  compares  the  experimental  and  simulated  dis¬ 
charge  curves  at  selected  cycles  and  the  model  parameters  are 
listed  in  Table  1. 

Fig.  3  shows  the  change  of  the  three  model  parameters  with 
cycling.  The  anode  and  the  cathode  become  consistently  less 
charged  with  cycling  (vo,p0s  increases  and  vo,neg  decreases).  It 
is  more  severe  for  the  anode  than  for  the  cathode.  The  volume 
fraction  of  the  cathode  £pos  decreases,  indicating  that  the  cathode 
loses  active  material  with  cycling. 

Fig.  4  shows  the  predicted  values  of  SOCs  for  the  anode 
(xneg)  and  the  cathode  (xpos)  at  the  end  of  discharge  (EOD)  for 
each  cycle.  A  two- stage  pattern  is  clearly  seen  in  Fig.  4  for  the 
electrode  SOCs.  The  cathode  becomes  less  intercalated  at  EOD 
in  the  first  stage,  then  more  intercalated  in  the  second  stage. 
Although  the  anode  SOC  changes  in  a  similar  pattern,  it  must  be 


Time  (s) 

Fig.  2.  Comparison  of  the  experimental  and  simulated  discharge  curves  at 
selected  cycles.  The  symbols  are  experiment  data  and  the  lines  are  the  simulation 
predictions.  The  discharge  current  was  1  mA. 


Table  1 


Model  parameters  used  in  the  study 


Parameters 

Anode 

Cathode 

Maximum  solid  phase  concentration 

3.056E— 2 

5.156E— 2 

(mol  cm  3),  cmax,/ 

Solid  phase  diffusion  coefficient 

IE— 11 

2E-10 

(cm2s  l),Dsj 

Rate  constant  (cm2  5  mol-0  5  s-1), 

1.07E-4 

6.36E-5 

Molecular  weight  (g  mol-1),  MW; 

72 

97.68 

Particle  size  (cm), 

5E-4 

5E-4 

Electrode  volume  (cm3),  V; 

4.1 

3.5 

Initial  SOC  of  the  electrode,  xo,;a 

0.823 

0.2925 

Volume  fraction  of  active  material,  £;b 

0.3749 

0.34 

Electrolyte  concentration  (mol  cm-3),  ce 

1.0E- 

-3 

a  Initial  SOCs  for  the  anode  and  the  cathode  would  change  during  cycling. 
The  value  for  the  first  cycle  is  shown  in  the  table. 

b  The  volume  fraction  of  active  material  for  the  cathode  decreases  with  cycling. 
The  value  is  shown  for  the  first  cycle. 


Fig.  3.  The  change  of  model  parameters  with  cycles  (obtained  from  fitting  the 
discharge  data  for  each  cycle).  Both  the  anode  and  the  cathode  become  less 
charged,  as  xo,neg  decreases  and  xo,pos  increases.  The  volume  fraction  of  the 
cathode  £pos  decreases,  indicating  that  the  cathode  loses  some  active  material 
with  cycling. 


Fig.  4.  The  simulated  SOCs  for  the  anode  and  the  cathode  at  the  end  of  discharge 
(EOD)  predicted  by  the  model  using  parameters  shown  in  Fig.  3.  A  two-stage 
behavior  is  clearly  seen  from  the  figure  for  the  electrode  SOCs. 
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Li1.xNiCo02 

0.0  0.2  0.4  0.6 


Fig.  5.  The  change  of  stoichiometric  windows  (or  SOCs)  for  the  electrodes  with 
cycling.  The  cross  marker  (X  )  indicates  the  movement  for  xo,pos  and  xo,neg  and 
the  circle  marker  (O'  )  indicates  the  movements  for  xpos  and  xneg  (SOC  at  the 
EOD).  Note  the  cathode  OCP  is  plotted  against  1  —  x  in  Lii-^NiCoCF.  Both 
electrodes  become  less  charged  with  cycling.  A  two-stage  pattern  is  observed 
for  the  electrode  SOC  at  the  EOD  (refer  to  Fig.  4  for  xpos  and  xneg). 

noted  that  the  anode  SOC  at  EOD  changes  on  a  much  smaller 
scale,  between  0.042  and  0.045,  as  compared  to  the  cathode 
SOC  at  EOD.  A  detailed  discussion  on  the  two-stage  pattern  is 
presented  below.  Movements  of  the  stoichiometric  windows  for 
the  electrodes  are  summarized  in  Fig.  5. 

3.2.  Model  prediction 

The  model  parameters  are  extrapolated  in  Fig.  6  (see  Table  2). 
The  assumption  is  that  the  capacity  fade  mechanisms  remain 
unchanged  during  the  period  of  our  prediction. 


Fig.  6.  The  extrapolation  of  the  model  parameters  shown  in  Fig.  3.  The  expres¬ 
sions  used  for  the  extrapolation  are  listed  in  Table  2. 


Table  2 


The  extrapolation  expressions  for  the  parameters  shown  in  Fig.  6 


Parameters 

Expression 

a 

b 

c 

-V),pos 

ax8  +  c 

3.656E— 5 

1.221 

0.2928 

•V),neg 

ax8  +  c 

— 4.262E— 3 

0.7507 

0.8267 

£pos 

aid*  +  c 

— 3.468E— 4 

1.053 

0.3406 

Cycles 

Fig.  7.  The  discharge  capacities  predicted  by  the  model  compared  to  the  one 
predicted  by  the  empirical  expression.  A  stage  of  accelerated  capacity  fade  is 
predicted  by  the  physics-based  model. 

The  predicted  discharge  capacity  as  a  function  of  cycle  num¬ 
ber  is  shown  in  Fig.  7  using  the  extrapolated  model  parameters. 
A  stage  of  accelerating  capacity  fade  is  predicted  after  about 
390  cycles.  Also  shown  in  Fig.  7  is  a  prediction  of  the  discharge 
capacities  obtained  by  using  an  empirical  equation  (see  Table  3). 
It  should  be  mentioned  that  researchers  [4]  have  observed  capac¬ 
ity  fade  patterns  shaped  in  a  similar  means  to  that  shown  in 
Fig.  7. 

Fig.  8  shows  the  simulated  electrode  SOCs  at  the  EOD  using 
the  extrapolated  model  parameters.  Please  note  that  the  scale  for 
Xneg  in  Fig.  8  is  different  from  the  one  used  in  Fig.  4.  A  third 
stage  is  observed  for  the  change  of  xpos  and  xneg  with  cycling. 
The  SOC  of  the  cathode  at  the  EOD  becomes  mostly  fixed  at 
0.98  in  this  stage,  indicating  that  the  cathode  is  close  to  being 
fully  intercalated.  The  SOC  of  the  anode,  however,  increases 
dramatically  (less  discharged)  at  the  EOD,  as  compared  to  the 
negligible  movement  in  the  first  two  stages.  Fig.  9  shows  the 
movement  of  the  SOCs  for  the  electrodes  in  all  the  three  stages. 

By  comparing  Figs.  6  and  7,  it  is  found  that  the  appear¬ 
ance  of  the  three  stages  for  the  electrode  SOCs  at  the  EOD 

Table  3 


Parameter  values  for  the  empirical  expression  used  in  Fig.  7 


Expression 

Ax8  +  Cx  +  D 

A 

-0.00971 

B 

0.5007 

C 

-8.03  x  10“3 

D 

0.9946 
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«  0.91 
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Cl 
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0.86 


0.81 

0  200  400  600 

Cycles 

Fig.  8.  The  simulated  SOCs  for  the  anode  and  the  cathode  at  the  EOD  using 
extrapolated  model  parameters  shown  in  Fig.  6.  The  symbols  represent  the  SOC 
data  plotted  in  Fig.  4.  Please  note  the  scale  for  xneg  is  different  from  the  one  in 
Fig.  4.  A  third  stage  shown  for  xpos  and  xneg  is  predicted  by  the  model. 
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corresponds  to  the  different  capacity  fade  stages  shown  in  our 
discharge  capacity  predictions.  Next,  we  will  analyze  the  capac¬ 
ity  fade  of  the  cell  based  on  the  information  obtained  using  the 
physics-based  model. 

3.3.  Capacity  fade  analysis 

Our  study  shows  that  there  could  be  multiple  stages  in  the 
capacity  fade  of  the  cell.  In  this  section,  the  scenarios  behind 
those  stages  for  cell  capacity  fade  are  presented. 

The  parameter  estimation  study  shows  that  the  volume  frac¬ 
tion  of  the  cathode  decreases  with  cycling  (£pos  in  Fig.  3), 
indicating  that  the  cathode  loses  some  of  the  active  material 
in  the  electrode.  In  addition,  it  is  well  known  that  there  are 
side  reactions  happening  on  the  anode  during  charge  process, 
which  would  consume  the  cyclable  lithium  ions  to  form  SEI 


Li^iCoOj 


0.0  0.2  0.4  0.6 


Fig.  9.  The  simulated  SOCs  for  the  electrodes  at  the  EOD  using  extrapolated 
model  parameters. 


film  on  the  anode.  The  major  capacity  fade  mechanisms  in  the 
cell  simulated  here  are  expected  to  be  the  loss  of  lithium  ions  to 
a  film  formation  reaction  and  the  loss  of  active  material  in  the 
cathode. 

The  loss  of  lithium  ions  to  a  film  formation  reaction  controls 
the  capacity  fade  in  the  first  stage,  because  the  film  formation 
reaction  rate  is  relatively  high  on  a  fresh  cell.  It  is  estimated  that 
during  the  first  charge/discharge  cycle,  about  18%  of  the  charge 
capacity  is  irreversibly  lost  to  the  film  formation  [  14] .  In  our  case, 
the  cathode  is  charged  from  Li(i.o)NiCoC>2  to  Li(o.3)NiCo02 
(Affpos  ~  0.3  in  Fig.  3)  on  the  first  charge.  Then  the  cathode  is  only 
intercalated  back  to  Li(o.89)NiCo02  at  the  end  of  first  discharge 
(■Tpos  ~  0.89  in  Fig.  4).  The  loss  ratio  is  (1— 0.89)/(l— 0.3)  =  0.16. 
That  is,  16%  of  the  charge  capacity  (1-0.3)  is  irreversibly  lost 
to  the  film  formation  for  the  first  cycle,  which  is  in  good  agree¬ 
ment  with  the  18%  estimation  in  Ref.  [14].  After  the  first  cycle, 
the  side  reaction  rate  gradually  decays,  because  the  SEI  film  on 
the  anode  partially  impedes  the  further  occurrence  of  the  side 
reaction.  This  is  supported  by  the  observation  that  the  capacity 
fade  rate  decreases  initially,  and  the  decrease  of  vo,neg  sl°ws  with 
cycling. 

In  the  first  stage  of  capacity  fade,  the  anode  is  the  limiting 
electrode.  Here  the  limiting  electrode  refers  to  the  one  which  is 
fully  discharged  at  the  EOD  and  causes  the  cell  to  reach  end  of 
discharge  voltage.  The  circle  marker  (O  )  in  Fig.  9  shows  clearly 
that  the  anode  is  fully  discharged  at  EOD.  The  cathode,  on  the 
other  hand,  is  not  fully  intercalated,  because  some  lithium  ions 
have  been  irreversibly  consumed  by  the  film  formation  reaction. 
As  the  loss  of  lithium  ions  becomes  larger,  the  cathode  is  con¬ 
tinuously  less  intercalated  at  the  EOD,  leading  to  a  decreasing 
Xpos  as  shown  in  Fig.  8. 

In  the  second  stage,  the  controlling  capacity  fade  mecha¬ 
nism  shifts  from  the  loss  of  lithium  ions  to  the  loss  of  active 
cathode  material.  As  mentioned  previously,  the  film  formation 
reaction  rate  decays  in  the  first  stage  after  the  SEI  film  forms 
on  the  anode,  leading  to  a  slower  loss  of  lithium  ions.  The  loss 
of  active  cathode  material  then  gradually  outpaces  the  loss  of 
lithium  ions,  which  causes  the  cathode  to  become  more  inter¬ 
calated  at  the  EOD  (vp0s  increases  in  Fig.  8).  Fig.  8  shows  this 
transition  point  as  a  minimum  in  the  vpos  versus  cycles,  where 
the  loss  of  cathode  material  starts  to  outpace  the  loss  of  lithium 
ions.  However,  the  anode  is  still  the  limiting  electrode  in  the  sec¬ 
ond  stage.  Fig.  8  shows  that  the  anode  is  mostly  discharged  to 
the  SOC  below  0.05  in  the  second  stage,  which  can  be  regarded 
as  fully  discharged  at  the  EOD.  Although  the  SOC  of  the  cath¬ 
ode  (vpos)  at  the  EOD  increases  with  the  cycling,  it  is  not  the 
cathode  that  causes  the  cell  to  reach  the  EODV  in  the  second 
stage. 

In  the  third  stage,  the  capacity  fade  mechanisms  remain  the 
same  as  in  previous  stages.  However,  the  limiting  electrode  has 
been  shifted  from  the  anode  to  the  cathode.  Fig.  8  shows  that 
the  cathode  is  nearly  fully  discharged  (to  LipoNiCo02)  at  the 
beginning  of  the  third  stage.  The  anode  becomes  increasingly 
less  discharged  at  the  EOD.  More  and  more  “cyclable”  lithium 
ions  are  left  inside  the  anode  and  cannot  be  intercalated  back 
into  the  cathode.  The  holding  capacity  of  the  cathode  decreases 
further  as  the  cathode  loses  more  active  material.  In  return,  more 
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“cyclable”  lithium  ions  are  kept  inside  the  anode  at  the  EOD, 
causing  accelerating  capacity  fade. 

4.  Conclusions 

A  physics-based  single  particle  model  is  used  to  study  the 
cycling  data  of  a  lithium  ion  cell.  Parameter  estimation  technique 
is  applied  to  extract  model  parameters  from  the  experiment  dis¬ 
charge  data.  The  study  on  the  experiment  data  shows  that  there 
are  two  stages  in  the  cell  capacity  fade.  Model  simulation  fur¬ 
ther  indicates  that  there  is  a  third  stage  of  capacity  fade  yet  to 
come  where  the  limiting  electrode  shifts  from  the  anode  to  the 
cathode. 

The  loss  of  lithium  ions  to  a  him  formation  reaction  (SEI 
formation)  and  the  loss  of  active  cathode  material  are  the  major 
capacity  fade  mechanisms.  In  the  first  stage,  the  loss  of  lithium 
ions  to  the  him  formation  reaction  is  the  controlling  mechanism. 
However,  the  side  reaction  rate  decays  after  the  SEI  him  forms 
on  the  anode.  The  cathode  becomes  less  intercalated  at  the  EOD 
as  cyclable  lithium  ions  are  irreversibly  consumed  in  the  him 
formation  reaction.  The  anode  is  the  limiting  electrode  in  the 
hrst  stage.  In  the  second  stage,  the  loss  of  active  cathode  material 
outpaces  the  loss  of  lithium  ions  because  of  the  decay  of  the  him 
formation  reaction  rate.  The  cathode  becomes  more  intercalated 
at  the  EOD.  The  anode  is  still  the  limiting  electrode.  In  the  third 
stage,  the  limiting  electrode  shifts  from  the  anode  to  the  cathode. 
Because  the  anode  could  not  be  fully  discharged,  “cyclable” 
lithium  ions  are  left  inside  the  anode,  causing  the  accelerating 
capacity  fade. 
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